The onset of magnetic order in fcc-Fe films on Cu(lOO) 
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3 ■ Abstract 

^ I On the basis of a first-principles electronic structure theory of finite tem- 

en I perature metallic magnetism in layered materials, we investigate the onset of 

L».^ I magnetic order in thin (2-8 layers) fcc-Fe films on Cu(lOO) substrates. The 

nature of this ordering is altered when the systems are capped with copper. 

C^ I Indeed we find an oscillatory dependence of the Curie temperatures as a func- 

^^ , tion of Cu-cap thickness, in excellent agreement with experimental data. The 

jh : 

Q , thermally induced spin-fluctuations are treated within a mean-fleld disordered 

O ' 

local moment (DLM) picture and give rise to layer-dependent 'local exchange 

r> ' splittings' in the electronic structure even in the paramagnetic phase. These 

features determine the magnetic intra- and interlayer interactions which are 

strongly influenced by the presence and extent of the Cu cap. 
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The Disordered Local Moment (DLM) picture |T[, when implemented by the KKR-CPA 
method 0, is a first-principles, that is to say material specific and yet parameter-free, mean 
field theory of metallic magnetism in the limit where good local moments form. It works well 
for bulk Fe and Co and their alloys but appears to be inadequate for Ni 0Q. Under these 
circumstances it is of fundamental interest to test its applicability by using it in different, 
but nevertheless similar, physical situations such as in a thin films whose thickness, d, can 
act as a tuneable, new thermodynamic variable. With this in mind we deployed the theory 
for Fe on and embedded in Cu. In this letter we report some of our most significant findings. 

For complexity in magnetic systems ultrathin Fe films on Cm(IOO) are unsurpassed [^. 
How the films grow, their structure and morphology are profoundly intertwined with their 
magnetic properties and a satisfactory description of these systems has become a benchmark 
for theories of thin film magnetism. After much extensive experimental work on thermally 
deposited films it is now clear that below a critical thickness of 10-12 monolayers (ML) the 
films take on the fee structure of the substrate while thicker films revert to the bcc structure 
of bulk iron. In the ultra thin regime when there are fewer than 4-5 ML a ferromagnetic (FM) 
phase is observed whilst the thicker films of 6-11 ML seem to be antiferromagnetic (AF) with 
a net moment across the film. Adding covering layers of Cu to the films has a marked effect 
upon their magnetic properties. A single Cu monolayer suppresses the magnetic ordering 
temperatures T^. whilst the T^'s of 2ML Cu-capped films are partly restored. In fact T^ 
oscillates weakly as further Cu layers are added 0. 

The experimental data on the uncapped films has been interpreted |^ in terms of mod- 
els in which the top two surface iron layers are coupled together ferromagnetically in an 
otherwise AF Fe film. This picture has been prompted by a series of ab-initio T = OK 



electronic structure total energy calculations [§-|T0|. Pajda et al. |TT[] have calculated ef- 
fective exchange interactions in a Cw-capped Fe monolayer on Cu(lOO) and obtained an 
oscillatory variation of Tc with cap thickness. Here we describe the full DLM picture of 
both electronic and magnetic structure pertinent to the paramagnetic states of Fe-films on 
C'u(lOO) of up to 8ML thickness at finite temperatures and investigate both the growth of 



magnetic correlations and onset of magnetic order. Moreover we are able to extract intra- 
and interlayer 'exchange' interactions. 

The DLM picture is based on assuming a separation between fast and slow degrees of 
freedom in the interacting many-electron system. For times, r, long in comparison with 
an electronic hopping tiraes,h/W (~ 10^^^ seconds), where W represents a relevant band- 
width, but short when compared with typical spin fluctuation times, the spin orientations 
of the electrons leaving an atomic site are sufficiently correlated with those arriving that 
the magnetisation integrated over a unit cell and averaged over r is non-zero. These are the 
'local moments' which can change their orientations, described by a set of unit vectors {cj} 
labelled by the site index i, on a time scale longer than r while their magnitudes fluctuate 
rapidly on the time scale r. This timescale demarcation has been invoked in much recent 
work |[T^-[T5| on bulk systems but has not been examined and tested in thin film systems. 



This is the issue we address in this letter. 

For an ab-initio implementation of this picture, standard Spin Density Functional the- 
ory for studying electrons in spin-polarised metals is adapted to describe the states of 
the system for each orientational configuration {cj}. At the heart of the theory is the 
generalised electronic Grand Potential fijcj} of the system so constrained. A long time 
average can be replaced by an ensemble average with respect to the Gibbsian measure 
P{ei} = Z^^ exp —i3fl{ei}, where the partition function Z = Hi I diiexp —i3Q{ei}. The 
thermodynamic free energy, which accounts for the entropy associated with the orientational 
fiuctuations as well as creation of electron-hole pairs, is given hj F = —ksT log Z. Evidently 
f2{ej} plays the role of a classical 'spin' (local moment) Hamiltonian, albeit a highly compli- 
cated one. By choosing a suitable reference single-site 'spin' Hamiltonian fio{ei} = J2i ^i{^i) 
and using the Feynman-Peierls' inequality |T6[ a mean field theory is constructed [QQ . This 



is a 'first principles' formulation of the Disordered Local Moment (DLM) picture which can 
be implemented by an adaptation of SCF-KKR-CPA method ideally suited for calculating 
the partial averages uji{ei). The approach can be further improved via the construct of a 
generalised Onsager cavity field [|]. 



When the above procedure was used to study bulk magnetic metals, it was found that the 
local electronic structure can possess a 'local exchange' splitting even in the paramagnetic 
state [0]. This means that an electron spin-polarised parallel to a local moment will have 
a different density of states to that polarised anti-parallel. When all orientations of the 
moments are averaged over in the paramagnetic state the electronic structure is inevitably 
unpolarised but consequences from the presence of local moments can still be identified. 
(Moreover the magnitude of this splitting is expected to vary sharply as a function of wave- 
vector and energy. At some points, if the 'bands' are flat in a region of wave- vector and energy 
space, the local exchange splitting will be roughly of the same size as the rigid exchange 
splitting of the majority and minority-spin bands of the ferromagnetic state whereas at other 
points, where the 'bands' have greater dispersion, the splitting will vanish entirely.) This 
local exchange-splitting is the cause of the establishment of a local moment. Photoemission 
r?! and inverse photoemission (IPES) experiments [0 have found these qualitative features 



in bcc Fe. In our generalisation of the DLM theory to layered systems we find that the 
layer dependence of similar features in the electronic structure of paramagnetic Fe grown 
on a Cm(IOO) substrate drives the onset of magnetic order and the form of the magnetic 
interactions between Fe layers. 

To probe the content of the theory consider the response of the paramagnetic DLM 
state to the application of a small, external spin-only magnetic field, {hi}, varying from 
site to site, layer to layer in both orientation and magnitude. The induced magnetisation 
is predominantly caused by the local moments changing their orientations to align with 
the field causing a change in the single-site probabilities, (5Pj(ej). In short the site by site 
paramagnetic spin susceptibility is the solution of the matrix equation, 

x^^ = f/^.% + fE4V^' (1) 

(2) 

where S^^ is the corresponding direct correlation function and /ij is the magnitude of the local 

(2) 

magnetic moment on the i-ih site. Expressions for S^j involving the electronic structure 
of the paramagnetic state and techniques for calculating them, using the KKR-CPA, for 



bulk systems have been given elsewhere 0. For layered systems with two-dimensional 
translational symmetry the magnitudes of the 'local moments', {yUj}, assume one value fip 
per site in a given layer labelled P but vary from layer to layer. Taking a 2D lattice Fourier 
transform over sites within each layer equation (|1|) can be rewritten as, 

^''^(qil) = Ex'^^^^' exph^qii ■ (R,, - R,)] 



V 



= f/^P^PQ + fE4%||)x"«(q||) (2) 

where, qy is a wave- vector in the 2D layer Brilliouin Zone. Once the S'p5'(q||)'s have been 
obtained from the SCF-KKR-CPA calculations and loaded into a n x n matrix, S''^^^(q||) 
(where n is the number of layers in the film), this set of equations (0) is solved by a simple 
matrix inversion i.e.x'^*^(q||) = SksTI — S'*-^''(q||) /ig (/ is a unit matrix). For films in 
which the intralayer exchange coupling is ferromagnetic, as for the Fe films on Cu, the Curie 
temperature Tc is specified by the condition || SfcsTcJ — S'*^^^(q|| = 0)|| = or in terms of the 
largest eigenvalue of S'*^^)(q|| = 0). Full technical details on how ^(^^(qii) is calculated for 
layered systems will be provided elsewhere. 

We use the spin-polarised screened KKR method [|I^ for layered systems adapted for 
the DLM picture. The lattice constants of the layers are taken to be the same as that of 
the substrate (6.83 a.u.) so that the effects of lattice strain on the electronic structure are 
neglected. For each n, the electronic structure of the DLM state is calculated self-consistently 
using 78 k\\ points in the irreducible part of the surface Brillouin zone. A buffer of three 
layers of the substrate as well as a buffer of at least three layers of vacuum is calculated 
self-consistently along with the potentials on each layer. The self-consistent layer-resolved 
potentials are then used to calculate the magnetic moments and the effective 'Weiss' field on 
each layer as described above. By introducing a small change in the average magnetisation 
on one particular layer and calculating the effective 'Weiss' field on each layer we determine 
'S'pQ (qii = 0), and hence the static paramagnetic spin susceptibility, x^'^(q|| = 0) together 
with the Curie temperatures Tc's for various thicknesses of both the Fe film and the Cu 
cap. 



Table | shows our results for the layer-dependent 'local' magnetic moments in uncapped 
Fe„/Cu(100) films together with the dependence of the Curie temperatures on Fe-film thick- 
ness. The moments on the topmost layer are always the largest at around 2.5/iB followed by 
the moments on the Fe layer closest to the substrate next at roughly 2.2^b- In the interior of 
the film the moments are reduced to 1.7 ^b close to what was found for fcc-Fe for this lattice 
spacing [^ . The magnetic ordering temperatures reduce monotonically to a constant value 



of 485 K for 6 MLs onwards. Table || shows the effective 'exchange' interactions, that is to 

(2) 

say the direct correlation function, Spg which lead to these transition temperatures. Exam- 
ples of a 7ML and a 3ML Fe layer systems are shown to bring out the dominant features. 
The intralayer couplings are ferromagnetic throughout. The layer at the top of film and 
the layer nearest the Cu substrate have the largest values owing to their reduced effective 
coordination. It is the interlayer couplings that show an important and interesting trend. 
The top two layers are strongly ferromagnetically (FM) coupled (positive values) whereas 
the coupling to nearest neighbor layers within the films is anti-ferromagnetic (AF). There is 
also significant coupling to subsequent layers which alternates between FM and AF. 

Adding a copper cap to these films alters these results profoundly. Figure |I| highlights 
this for 3 and 7 monolayer Fe/Cu(100) films and shows the magnetic ordering temperature 
Tc to have an oscillatory dependence on cap thickness. A single Cu monolayer suppresses 
Tc by some 50K with the deficit being restored by a second layer. A third layer reduces T^ 
once again and further layers cause minor oscillations about this lower value. This behavior 
is in excellent agreement with experimental data on the same systems by Vollmer et al. 
(see their figure 6). In figure ^ the dominant changes to intra- and interlayer magnetic 
interactions are shown. These occur in the top two Fe layers nearest the cap. A single Cu 
ML cap switches the coupling between these layers from FM to AF whilst further Cu layers 
strengthen the coupling within the topmost Fe layer. 

Calculations for bulk fcc-Fe have shown how the tendency to magnetic order changes as 
the lattice spacing is varied pT|,pU|. On an expanded lattice, a ferromagnetically ordered 



state is stable at low temperatures and the coupling between 'local moments' in the param- 

6 



agnetic DLM state mirrors this aspect. As the lattice is contracted the magnetic ordering 
tendencies become anti-ferromagnetic. Capping fee Fe films with Cu causes a similar effect. 
Figure ^ shows that the electronic density of states (DOS) of an uncapped film is 'exchange 
split' over a large energy range. These states are nearly fully occupied and FM coupling be- 
tween neighboring sites results. The reduced coordination of sites within this layer produce 
an effect akin to that of increasing the lattice spacing of bulk fee Fe. There is hybridisation 
of these Fe states with quantum well states in the Cu cap formed between the vacuum and 
the Fe film. This causes the exchange splitting to collapse in the region of this hybridisation 
just as in fee Fe for smaller lattice spacings. The remaining exchange-split states are roughly 
half-filled which leads an AF coupling along the (100) direction as in bulk fee Fe. Figure |^ 
shows these effects in the DOS in the top layer of a Fe7/Cu(100) system both uncapped and 
capped by a single Cu ML. The narrow band of states in the Cu layer are also shown. 

In conclusion we note that the dramatic variations of interaction between local moments 
at different distances from the various surfaces and interfaces are quintessentially itinerant 
effects. As the much studied variation of the local moments near such planar defects show, 
they are due to the motion of the electrons not limited to an atom but confined by the 
geometry of the sample. In this respect they are manifestations of the same physics as is 
at work in the oscillatory magnetic coupling in metallic multilayers |2^. Consequently, the 



detailed quantitative account of the complex experimental data by our calculations can be 
taken as significant new evidence that the first-principles (KKR-CPA) implementation of 
the DLM picture correctly captures the essential physics of magnetic order, due to mobile 
electrons, in the limit where local moments form. 
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TABLES 
TABLE L Tocal' magnetic moments in fis on different layers and the Curie temperatures of 

uncapped Fe„/Cu(100) systems. The Fe Layer Li is adjacent to the substrate and layer L„ is the 

top layer. 





Tocal' moments {^b) 


Tc{K) 


n 


Li 


L2 


Ls 


L4 L5 


Le 


Li Lg 




2 


2.24 


2.48 










681 


3 


2.21 


1.70 


2.53 








532 


4 


2.21 


1.67 


1.74 


2.53 






495 


5 


2.21 


1.67 


1.71 


1.74 2.53 






492 


6 


2.21 


1.67 


1.71 


1.71 1.74 


2.53 




486 


7 


2.21 


1.67 


1.71 


1.72 1.71 


1.74 


2.53 


485 


8 


2.21 


1.67 


1.71 


1.72 1.71 


1.71 


1.74 2.53 


485 



TABLE IL Intra- and inter-layer effective 'exchange' interactions in meV between layers in an 
uncapped Fe7/Cu(100) system with the values for Fe3/Cu(100) shown for comparison in brackets. 
The Fe Layer Li is adjacent to the substrate and layer L„ is the top layer. 



n 


Li 


L2 


L3 


L4 


L5 


Le 


Ly 


1 


108.5 (108.4) 


-7.5 


8.5 


-2.1 


-0.8 


-0.7 (-10.9) 


-1.0 (25.1) 


2 


-7.5 


19.6 


-7.5 


3.2 


-2.2 


-0.6 


-0.7 


3 


8.5 


-7.5 


17.5 


-7.2 


2.1 


-1.7 


-0.7 


4 


-2.1 


3.2 


-7.2 


18.8 


-8.1 


2.1 


-2.1 


5 


-0.8 


-2.2 


2.1 


-8.1 


19.9 


-10.0 


9.8 


6 


-0.7(-10.9) 


-0.6 


-1.7 


2.1 


-10.0 


33.2(29.9) 


50.0 (53.8) 


7 


-1.0(25.1) 


-0.7 


-0.7 


-2.1 


9.8 


50.0 (53.8) 


97.7 (96.7) 
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FIG. 1. The magnetic ordering temperatures for 3 and 7 layers of Fe on Cu{WO) as a function 
of number of copper capping layers. 
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FIG. 2. The magnetic interactions in and between the Fe layers nearest the surface or copper 
cap in a system of 7 Fe layers on (7u(100). The values are given in meV and ^L^j-Lj^ denotes the 
interactions within the top Fe layer (nearest the cap), ^L^-Lq for the second-most top layer and 
^Lq-Lj^ the coupling between these two iron layers. 
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FIG. 3. The density of states (DOS) in the topmost Fe layer of a Fe7Cu(100) system. The full 
line shows the uncapped system whilst the dashed line shows the DOS when a single Cu monolayer 
caps the system. The dotted line shows the DOS for this Cu layer. The upper half of the figure 
shows the DOS for an electron spin-polarised parallel to the local moment on a site whereas in the 
lower half the DOS of an electron polarised anti-parallel is shown. 
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